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Abstract. Very weakly interacting slim particles (WISPs), such as axion-like particles (ALPs) 
or hidden photons (HPs), may be non-thermally produced via the misalignment mechanism 
in the early universe and survive as a cold dark matter population until today. We find that, 
both for ALPs and HPs whose dominant interactions with the standard model arise from 
couplings to photons, a huge region in the parameter spaces spanned by photon coupling 
and ALP or HP mass can give rise to the observed cold dark matter. Remarkably, a large 
region of this parameter space coincides with that predicted in well motivated models of 
fundamental physics. A wide range of experimental searches - exploiting haloscopes (direct 
dark matter searches exploiting microwave cavities), helioscopes (searches for solar ALPs or 
HPs), or light-shining-through-a-wall techniques - can probe large parts of this parameter 
space in the foreseeable future. 
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1 Introduction 

The two most relevant features of dark matter (DM) particles are their feeble interactions 
with standard model particles and their cosmological stability. In addition DM is required 
to be sufficiently cold in order to allow for efficient structure formation. 

A realization of all these features are weakly interacting massive particles (WIMPs). 
Thermally produced in the early Universe their large, of order TeV scale mass ensures that 
by now they are non-relativistic and therefore sufficiently cold. Their interactions are small 
due to the large mass of the mediator particles (such as W or Z bosons) which makes the 
interaction very short ranged. Despite this, their large mass and the correspondingly large 
phase space is at odds with the required stability on cosmological time-scales. In order to 
ensure this stability one is therefore forced to introduce symmetries that conserve the number 
of these DM particles. However, motivating these symmetries on theoretical grounds is non- 
trivial. Global symmetries may be broken in quantum gravity whereas local symmetries lead 
to additional interactions which may cause conflicts with the required weakness of the DM 
interactions. Nevertheless good candidates exist. Two of the most famous examples are: 
the lightest supersymmetric particle in supersymmetric models with R-parity (see, e.g. [1] 
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for a review), or the lightest Kaluza-Klein modes in models with conserved parity in extra 
dimensions (see, e.g. [2] for a review), both models very appealing because of their connection 
with more fundamental theories of space and time, and their discovery potential at the LHC. 

Although it is way too early to make a final judgement it is nevertheless noteworthy 
that the initial searches at LHC have not given any indication of the existence of WIMPs. 
Because of this and the mentioned theoretical issues it is interesting and timely to consider 
alternative ways to realize the essential features of DM. 

Sufficient stability of the DM particles can also be achieved by combining the weakness 
of their interactions with a sufficiently small mass. The latter drastically reduces the phase 
space (and the number/type of possible decay products) thereby increasing the lifetime. This 
is the road we want to pursue in this paper: we will concentrate on very weakly interacting 
slim particles (WISPs) as DM candidates. 

Yet, a thermally produced light DM candidate can run foul of the required coldness 
of DM and can prevent structure formation. More precisely, the free-streaming length (the 
distance a DM particle can travel before getting trapped in a potential well) would increase 
with decreasing mass, and therefore at some point these DM particles would be inconsis- 
tent with the existence of dwarf galaxies, galaxies, clusters, superclusters and so on. This 
argument can be used for instance to rule out standard neutrinos as DM. This reasoning is 
extremely powerful in light of the increasingly precise cosmological data at our disposal and 
even subdominant components of thermally produced light DM can be ruled out, the case of 
eV mass axions being a prime example (see, e.g. [3]). 

However, there are non-thermal means for producing sufficiently cold dark matter made 
of light particles. One of the most interesting is the misalignment mechanism, discussed 
mostly in the case of the QCD axion [4-6] or (recently) string axions [7-10] and the central 
topic of this paper. Very recently, this mechanism has also been proposed to produce cold 
dark matter (CDM) out of hidden photons (HPs) [11]. 

In this paper, we shall revisit the misalignment mechanism for both cases: we treat 
first axion-like particles (ALPs - which may arise as pseudo-Nambu-Goldstone bosons in 
field theory and appear generically in all string compactifications) and then hidden photons. 
Our conclusions turn out to be extremely motivating. Once produced, a population of very 
light cold dark matter particles is extremely difficult to reabsorb by the primordial plasma. 
Therefore, we find that in both cases, ALPs and HPs, a huge region in parameter space 
spanned by their masses and their couplings to standard model particles can give rise to the 
observed dark matter. These topics have already been discussed in some detail in [11, 12]. 
The novelty in this work is that we shall provide new constraints and expose interesting 
regions of parameter space relevant for direct and indirect searches. 

The outline of this paper goes as follows: in Sect. 2 we review the essentials of the 
misalignment mechanism. In Sects. 3 and 4 we elaborate on two particular cases, axion-like 
particles and hidden photons, respectively. We discuss the cosmological constraints, also 
noting some misconceptions in the results of [11]. In Sect. 5 we discuss the direct detection 
of ALP and HP CDM in microwave cavity experiments. We conclude in Sect. 6. 

2 Essentials of the misalignment mechanism 

The misalignment mechanism relies on assuming that fields in the early universe have a 
random initial state (such as as one would expect, for example, to arise from quantum 
fluctuations during inflation) which is fixed by the expansion of the universe; fields with mass 
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m evolve on timescales t ~ m~ . After such a timescale, the fields respond by attempting 
to minimize their potential, and consequently oscillate around the minimum. If there is no 
significant damping via decays, these oscillations can behave as a cold dark matter fluid since 
their energy density is diluted by the expansion of the universe as p oc a~^, where a is the 
universe scale factor. 

In order to be more quantitative, let us revisit this mechanism for the simple case of a 
real scalar field with Lagrangian, 

C=^d^cPd^cP-^my + Ci, (2.1) 

where Cj encodes interactions of the scalar field with itself and the rest of particles in the 
primordial bath. We assume that the universe underwent a period of inflation at a value of 
the Hubble expansion parameter H = dloga/dt larger than the scalar mass. After inflation 
the field shall be approximately spatially uniform and the initial state is characterized by a 
single initial value, After inflation a period of reheating occurs, followed by a period of 
radiation dominated expansion. The equation of motion for </> in the expanding Universe is 

4> + 3H(j) + ml(l) = 0. (2.2) 

In general, the mass receives thermal corrections from Cj which might be crucial, thus = 
m^{t) should be understood. 

The solution of this equation can be separated into two regimes. In a first epoch, 
'iH S> m^, so (j) is an overdamped oscillator and gets frozen, = 0. At a later time, ti, 
characterized by 'iH{ti) = m(p{ti) = mi, the damping becomes undercritical and the field 
can roll down the potential and start to oscillate. During this epoch, the mass term is the 
leading scale in the equation and the solution can be found in the WKB approximation, 

"^"^^'i^) ^^^{l^fdt), (2-3) 

where (pi ~ (pi since up to ti the evolution is frozen. Note that in obtaining this solution we 
have not assumed a particular form for H but just its definition, and so it is valid for the 
radiation, matter, and vacuum energy dominated phases of the universe and their transitions. 

The solution corresponds to fast oscillations with a slow amplitude decay. Let us call 
this amplitude A{t) = (/)i(miof /m^o^)^/^ and the phase a{t) = m^{t)dt. The energy 
density of the scalar field is 

P<t> = + \^\<i>^ = l^lA^ + - > (2-4) 

where the dots stand for terms involving derivatives of A, which by assumption are much 
smaller than m,^ (m^ ^ in this regime). Let us also consider the pressure, 

P(p = 2*^^ ~ 2"^^*^^ ~ "2"^^"^^ ^'^^ ^^'^•^ ~ AArUfp sin (2a) + A'^ cos^ (a) . (2.5) 

At times t ^ ti^ the oscillations in the pressure occur at time scales much much faster 

than the cosmological evolution. We can therefore average over these oscillations, giving 

(p^) = (^2^os2(a)) = ^^2_ (2.6) 

^When the field just starts to oscillate the averaging employed in the following is not a good approximation 
and the equation of state is a non-trivial and strongly time dependent function. Depending on when the 
transition occurs this may have interesting cosmological effects on, e.g., structure formation. 
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As already mentioned A <^ m^A. Thus, at leading order in A/ {Am), the equation of state 
is just 

w = {p)/{p) ~ 0, (2.7) 

which is exactly that of non-relativistic matter. 

It follows from (2.3) that the energy density in a comoving volume, pa^, is not conserved 
if the scalar mass changes in time. The quantity 

N = pa^/m^ = ^mial(j)l (2.8) 

is however constant even in this case, and can be interpreted as a comoving number of non- 
relativistic quanta of mass m^jj. Here, we only need the conservation of to compute the 
energy density today as 

P<l>{to) = ~ -momicpi — , (2.9) 

where quantities with a 0-subscript are evaluated at present time. 

More physics insight is gained using temperatures instead of times and scale factors. 
First, we use the conservation of comoving entropy S = sa^ = 2Tr g^: s {T)T^ a"^ / 45 to write 
(ai/a)^ = g^s{T)T'^ / g^fs{Ti)Tf . Then we use the expression for the Hubble constant in the 
radiation dominated era H = l.m^Jg^{T)T'^/m^l and the definition of Ti, ZH{Ti) = 772-1 to 
express Ti in terms of mi and the Planck mass mpi = 1.22 x 10^^ GeV. The functions 5* and 
g-tS are the effective numbers of energy and entropy degrees of freedom defined in [13]. The 
dark matter density today, (2.9), can then be expressed as 

where F{Ti) = (g*(Ti)/3.36)4(g^5'(Ti)/3.91)~^ is a smooth function ranging from 1 to ~ 0.3 
in the interval Ti G (To,200GeV). The abundance is most sensitive to the initial amplitude 
of the oscillations, oc 4>\, and to a lesser degree to today's mass niQ. The factor oc 1/y/mi 
reflects the damping of the oscillations in the expanding universe: the later the oscillations 
start, i.e. the smaller Ti and therefore Hi and mi, the less damped they are for a given mo- 
If we compare the above estimate with the DM density measured by WMAP and other 
large scale structure probes [14], 

PCDM = 1.17(6)^, (2.11) 

cm'^ 

it is clear that we need very large values of (pi to account for all the dark matter. However, 
a relatively small (j)i could be compensated by a small mi ^ mo- 

If we want the condensate to mimic the behaviour of standard cold dark matter we 
should ensure that, at latest at matter-radiation equality, at a temperature Tcq ~ 1.3 eV, the 
mass attains its current value mo and therefore the DM density starts to scale truly as 1/a^. 
In particular, at this point the field should already have started to oscillate. This corresponds 
to a lower limit^ on mi, mi > 3if(req) = 1.8 x lO"'^^ eV, which implies an upper bound on 

keV mo / 4>i 
OT? ^ eV V 53TeV , 

2 See also [151. 
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In other words if we want these particles to be the DM, we need that (mo/eV)((/)i/53 TeV)^ > 
1, giving us a constraint on the required initial field value as a function of the mass today. 

To conclude this section let us note that dark matter generated by the misalignment 
mechanism may have interesting properties beyond those of cold dark matter. At the time 
of their production, particles from the misalignment mechanism are semi-relativistic. Their 
momenta are of the order of the Hubble constant p ~ Hi <^ Ti ; accordingly we have (outside 
of gravitational wells) a velocity distribution with a very narrow width of roughly, 

5v{t) ~ ^ f « 1. (2.13) 
mi \aoJ 

Combined with the high number density of particles, n^^g = -^/^o ~ /'cdm/^t-O) this narrow 
distribution typically leads to very high occupation numbers for each quantum state. 

A/occupation ~ 777 — ~ lU , (^•14) 

47r/3 rnQdv-^ \mo J \mo J 

where we have used oo/ai ~ Ti/Tq ~ ^/m^mpi/TQ. If the interactions are strong enough to 
achieve thermalisation, as argued in Refs. [16, 17] for the case of axions, this high occupa- 
tion number can lead to the formation of a Bose-Einstein condensate. This could lead to 
interesting properties which may also lead to interesting signatures in cosmological observa- 
tions [16-20]. Although we will not investigate this intriguing possibility here, we note that 
these features could also be realized for the more general light DM particles discussed in this 
paper. 

In the following we discuss two particularly interesting possibilities for DM from the 
misalignment mechanism, axion-like particles and hidden photons. 



3 Axion-like particles 

In this section we will focus on a specific type of WISP, namely axion-like particles (ALPs). 
By this we shall mean particles with only derivative couplings to matter, and in particular 
an interaction with photons given by 

C D -l^^F^^F/^^ (3.1) 

where (p is the ALP and g is a dimensionful coupling. The chief examples of ALPs are pseudo 
Nambu-Goldstone bosons (pNGBs) and string "axions" which can be treated together. For 
concreteness we will focus in this paper on particles of these types. 

The cosmology of the ALP condensate depends on the type of interaction generating 
its mass and in particular how this mass changes through the evolution of the universe. In 
the following we will go through a variety of possibilities for the DM formation as well as its 
cosmological viability for different scenarios. The regions which allow for viable ALP DM 
are then assembled in Fig. 1 in the — g plane. 

3.1 Axion-like particles from pNGBs and string theory 

When a continuous global symmetry is spontaneous broken, massless particles appear in the 
low energy theory: Nambu-Goldstone bosons (NGBs). They appear in the Lagrangian as 
phases of the high energy degrees of freedom. Since phases are dimensionless the canonically 
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normalised theory at low energies always involves the combination (j)/ f^, where (j) is the 
NGB field and is a scale close to the spontaneous symmetry breaking (SSB) scale. The 
range for (j)/ fcf, is (— vr, vr) and therefore the natural values for are ~ String axions on 
the other hand appear in all compactifications. They share these properties (having a shift 
symmetry and being periodic) but with the natural size of being the string scale (in type 
II compactifications this can be somewhat modified by a factor of the typical length scale of 
the compactification) . 

Indeed, all of the global symmetries in the standard model are broken^. Furthermore 
the black hole no hair theorem and what we know about quantum gravity tell us that this 
should ultimately also occur to any additional global symmetries. Hence we should have 
pseudo Nambu-Goldstone Bosons (pNGBs) instead of NGBs. They then have a mass, and 
can be a dark matter candidate. 

There are many possibilities for breaking the shift symmetry, explicitly or spontaneously, 
perturbatively or non-perturbatively; for stringy axions, the shift symmetry is exact to all 
orders in perturbation theory and is only broken non-perturbatively, for instance from a 
non-abelian anomaly, gaugino condensation or stringy instantons. The ALP potential can 
typically be parametrized as 

y(<A)=m2/2(^l-cos-^). (3.2) 

The mass of the ALP is in general unrelated to the QCD axion mass and in particular will be 
independent of the temperature unless generated by a sector that is thermalised. The ALP 
will satisfy the equation of motion (2.2) as long as (/"//^ is small, i.e. few oscillations after 
ti . The inaccuracy of the quadratic approximation can be cured by an additional correction 
factor to (2.10). This is normally an 0(1) factor except if we fine tune the initial condition 
to <p = 7r/0. 

The dimensionful coupling parameter g in (3.1) can be parametrised as 

In the simplest case Af is an integer, but this is not true in general when the ALP mixes, 
either kinetically or via symmetry-breaking effects with other ALPs or with pseudoscalar 
mesons. 

We can then represent the allowed regions of ALP dark matter in the m^-g plane by 
using 

(3.) 

with 9i = 1 01 1 //,/,, the initial misalignment angle whose range is restricted to values be- 
tween — vr and vr. The model dependent factor M will from now on be taken to be unity, 
but the reader should keep in mind that in principle it can be very different in particular 
constructions. 

While we shall adopt a phenomenological approach, showing the allowed region in the 
"^(/) ~ 9 plane, as mentioned above we are motivated by both field theoretical and stringy 
axion models, and so we may ask what the preferred region of the parameter space is. In 
both cases, the ALP could be related to the generation of right-handed neutrino masses, 

^Assuming that neutrinos are Majorana fermions, otherwise B — L is an exception. 
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and hence have a decay constant at an intermediate scale, i.e. within the cosmological 
axion window; alternatively it could be associated with a GUT, and have a decay constant 
at the corresponding scale. For string axions, the coupling to photons will be related via 
a loop factor to either the string scale or the Planck scale (or it could be even weaker - 
so an ALP with a large coupling would restrict the string scale to be low). For string 
and field theoretical models the most interesting values are therefore g ~ (10~^^, 10^^^) 
GeV~^, ~ 10~^^ GeV~^ and ~ 10^^^ GeV~^ corresponding to intermediate, GUT or Planck 
scales. On the other hand, we have no preconceptions regarding preferred ALP masses; as 
mentioned above the mass is expected to be generated by non-perturbative effects, and 
thus suppressed exponentially relative to some fundamental scale (in IIB string theory it 
is roughly the gravitino mass which multiplies the exponential factor) so for example with 
many stringy axions (an "axiverse" [7-10]) we could easily populate the allowed parameter 
region. In the following we will therefore treat both g and rUfp as theoretically unconstrained 
parameters. 

3.2 ALP dark matter from the misalignment mechanism 

The value of (pi that determines the DM abundance depends on the behaviour of the ALP 
field during inflation. For a pseudo-Nambu-Goldstone boson, the spontaneous symmetry 
breaking (SSB) could take place before or after inflation: the pNGB effectively exists only 
after SSB and it is during the associated phase transition that its initial values are set. For 
a string axion, provided the inflationary scale is below the string scale (or, equivalently, the 
decay constant) we should have control over the field theory, and so it will behave like a 
pNGB with symmetry broken before inflation. Assuming its mass to be much smaller than 
the Hubble scale at the time of SSB, Hssb, the ALP field will take random values in different 
causally disconnected regions of the universe. The initial size of these domains cannot be 
larger than 

1 mpi 
-Hssb f^y/g*{U) 

For a string axion, or for a pNGB whose associated SSB happens before inflation, 
the initial domains are stretched over a size larger than the current size of the universe. 
Consequently the initial field value is the same for every point within our current horizon. The 
current DM density then depends on this initial field value, leaving an additional parameter 
to tune the DM density. 

On the other hand, if the SSB happens after inflation, the DM density has inhomo- 
geneities of order 0(1) at length scales < Lj,dom- Non-linear effects, due to the attractive 
self-interaction caused by higher order terms in the expansion of the potential (3.2), drive the 
overabundances to form peculiar DM clumps that are called miniclusters [21-24]. These mini- 
clusters act like seeds enhancing the successive gravitational clumping that leads to structure 
formation. The minicluster mass is set by the dark matter mass inside the Hubble horizon 
dn = H^^ when the self-interaction freezes-out, i.e. Mmc ^ p^{T\)dH{T\)^ for the freeze-out 
temperature T\. Long-range interactions will be exponentially suppressed at distances longer 
than 1 /m^ so we can expect Tx to be of the order Ti , with at most a logarithmic dependence 
on other parameters. This is indeed the case for QCD axions, for which the minicluster- 
ing is quenched soon after the QCD phase transition that turns on the potential (3.2) [16] 
giving Mine ~ IO^^'^Mq, where Mq = 1.116 x 10^^ GeV is the solar mass, and a radius 
Rmc ~ 10^^ cm [25]. In the case of ALPs, can be larger if the mass is lighter. The 

authors of [24] pointed out that the present data on the CDM power spectrum constrain 
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Mrac ^ 4 X IO^Mq which translates into a lower bound in temperature > 2 x 10 ^ GeV 
and in the ALP mass > H{T = 2 x 10"^ GeV) ~ IQ-^o eV. 

If some of these miniclusters survive the tidal disruption during structure formation 
they should be observable in forthcoming lensing experiments [24, 25]. In any case, at larger 
scales, the DM density averages to a constant value corresponding to ~ 7r^/|/3 bearing 
the mentioned 0(1) correction due to the non- harmonic behaviour of large initial phases. 

During the spontaneous symmetry breaking of a global symmetry topological defects 
such as cosmic strings [26] and domain walls can be formed. Strings have a thickness (5 ~ 1//^ 
and typical sizes of the order of the horizon, L ~ t. As strings enter into the horizon they can 
rapidly reconnect, form loops and decay into pNGBs. If the SSB happens after inflation, we 
have to consider also their contribution to the energy budget of the universe. Axions resulting 
from string decay are known to contribute significantly to their cold DM density, but the exact 
amount is subject to a long-standing controversy [13]. The debate is focused around the axion 
emission spectrum. Some authors argue that the decay proceeds very fast, with an emission 
spectrum 1/k with high and low energy cutoff of order respectively 1/6 and 1/L. In this case 
the contribution to is similar to that from the misalignment mechanism [27, 28]. Others 
put forward that the string decays happen after many oscillations, with a radiation spectrum 
peaked around 27r/L, which enhances the contribution to cold DM by a multiplicative factor 
of log(L/5) ~ log{fa/ma) ~ 0(60) [29-35]. Once the axion potential builds up at the QCD 
phase transition, domain walls build up. If the axion potential has only one minimum these 
domain walls can still efficiently decay into axions, leading to a third axion population which 
is thought not to be significant in the fast decay scenario, while a recent analysis has shown 
that it can be the main contribution to the axion energy density in the many oscillations 
one [36]. If different exactly degenerate vacua exist the domain walls are persistent and 
can very easily run in conflict with observations. Therefore one assumes a small explicit 
breaking of the Peccei-Quinn symmetry, which breaks the degeneracy and allows domain 
walls to decay. It is possible although somehow fine-tuned to do so without compromising 
the solution to the strong-CP problem. For a recent review on axion cosmology see e.g. [37]. 

We expect the same type of behavior for ALPs with similar characteristics than axions, 
i.e. ALPs whose mass is generated at a late phase transition due to a hidden sector which 
becomes strongly interacting. In this case we should keep in mind the controversy of the string 
decay spectrum and assume an uncertainty of order log{ f^/in^) in the DM abundance. The 
domain wall problem can in principle be solved by strong enough explicit breaking, and their 
contribution to the DM appears to be subdominant as well. In models where the cosmological 
evolution of the ALP mass is different, the situation can differ from the above. These models 
have to be studied case by case, which is beyond the scope of this paper. 

3.3 Sufficient production of dark matter 

As we have seen in Sect. 2 a general constraint arises from the fact that we get a sufficient 
amount of DM but at the same time the mass at matter radiation equality has to be greater 
than the Hubble constant. This is the bound Eq. (2.12). For pNGB ALPs we however also 
have that the field value itself cannot be larger than irf^ which itself is connected to the 
coupling to photons. Combining these two restrictions gives us a way to constrain the viable 
regions. The light red region in Fig. 1 labelled "mi > 3H{Teq)" corresponds to this general 
bound with (pi < and J\f=l. 

Any ALP model should satisfy this bound for its zero mode to behave as DM before 
matter-radiation equality. Realistic models attempting to saturate this bound will have 
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Figure 1. Parameter space for axions (shaded band labelled "Axion models") and axion-like par- 
ticles. The regions where they could form DM are displayed in different shades of red (for details 
see text). The lines representing DM regions are uncertain through a model-dependent multiplicative 
factor, A^, which we have set equal to 1 here. The DM regions move towards larger couplings g, pro- 
portional to this factor. The exclusion regions labelled "ALPS", "CAST-f Sumico" and "HB" arise 
from experiments and astrophysical observations that do not require ALP dark matter (for a review, 
see [38]). The remaining constraints are based on ALPs being DM and are described in the text. 

problems either fitting the cosmic microwave background (CMB) data or with the WKB 
approximation we have used. In this sense this bound is very conservative. Importantly, for 
A/" ^ 1, it seems to exclude the possibility of providing DM from the type of ultralight ALP 
field that has been invoked to explain the puzzlingly small opacity of the universe for ~ TeV 
gamma rays (see Ref. [39] and references therein) in terms of photon -f-T- ALP conversions in 
astrophysical magnetic fields, requiring^ g ~ lO^^^GeV^^ and < 1 neV, see [40-43]. To 
allow an ALP to explain these observations and simultaneously to be dark matter requires 

*The required coupling is determined by the extragalactic background light and the size of the astrophysical 
magnetics fields, therefore plagued by sizeable uncertainties. 
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A/" > 10 which is still conceivable, or a fine tuning of 9i. 

Let us now turn to the stronger constraints that can be obtained for specific 
time/temperature dependencies of the ALP mass. The simplest realization of an ALP model 
has constant throughout the universe expansion. In this case, we can infer the DM yield 
from Eq. (2.10) using mi = uiq. The region of ALP DM in this case is depicted in pink 
in Fig. 1 and labelled "Standard ALP CDM (mi = mo)". We have assumed J\f = 1 and 
used the a-priori unknown value of 9i to tune the right DM abundance. The upper bound 
on g reflects the fact that 6i cannot be larger than vr, and thus assumes ~ vr. Moving 
to lower values of g requires inflation happening after SSB in order to have a homogeneous 
small value of 6i which is increasingly fine-tuned to zero. In this sense the values closest to 
the boundary, corresponding to the largest values of the photon coupling, can be considered 
the most natural values. 

Quantum fluctuations generated in the ALP field during inflation will produce unavoid- 
able inhomogeneities in 6i of order 66i ~ Hj / {2Trf(f,). This precludes fine-tuning of the 
universe average of 4>i below Hj /27r, and sets a minimum DM abundance for a specified 
value of Hj . 

Since the ALP is effectively massless during inflation, these inhomogeneities correspond 
to isocurvature perturbations of the gravitational potential; this has been discussed exten- 
sively in the literature in the context of axions (see e.g. [44]) and for many string axions, 
see [8]. The WMAP7 observations of the primordial density fluctuations set very stringent 
constraints on isocurvature perturbations from which one can obtain an upper bound on Hj 
assuming a given ALP DM defined essentially by f,f,. Assuming that the decay constant 
does not change during inflation (it certainly should not for string axions, for example) the 
constraint is [14] 

where (15*^1) is the isocurvature power spectrum, and the adiabatic one (generated by 

the inflation or other fields). We can approximate (I*? |) ~ where Hj is the Hubble 

constant during inflation. At the pivot scale ko = 0.002 MPc"^ WMAP finds {\R'^\) = 
2.42 X 10~^ so we have a bound 

Hi <4.5 X W~^(Pi. (3.7) 

In principle we can constrain the mass of the ALP at the time when it starts oscillating by 
requiring 3Hj > 3H{Ti) = mi, 



h >7 X lO^mi, 



.>...-.^.(HIM^)^Q'(^)\ (3.S) 

where in the last equation we used (2.10) to get an expression for mi. Therefore, at the 
boundary of sufficient dark matter production when cpi = vr/,^, and taking g^ = a/27r/0 this 
translates into 

1 >0.08^^ X ^"'"-'-'y Qy ('"'" Q"^-')'. (3.9) 
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Clearly this is a rather weak constraint, and we are therefore allowed many orders of mag- 
nitude between the inverse of the inflationary Hubble constant and the time when the ALP 
oscillations begin. Conversely, the reheating temperature Trh is bounded by T^jj < HjMp; 
with additional assumptions this can constrain 0i: for example, requiring leptogenesis with 
Trh ~ 10^ GeV would bound 0i > 10^ GeV. 

We do not have strong arguments against values of Trh larger than a few MeV [45] , only 
appealing theoretical prejudices at most, so we shall leave this question aside. In any case, 
the most interesting region from observational purposes corresponds to the largest values of 
where the initial angle is not fine-tuned. 

This simple ALP model predicts quite weakly coupled ALP CDM. Models in which 
mi ^ mo can provide larger DM abundances with smaller values of and therefore stronger 
interactions. 

An interesting example of this case arises if the ALP acquires a mass from coupling to 
some (hidden) non-abelian group. For a pNGB, this would mean that the associated global 
symmetry is anomalous under some the non-abelian group, just as the r]' or the hypothetical 
axion acquire their mass via QCD instantons. For our ALP we need in principle another 
unbroken SU(N) group, which condenses at a scale A. Then we can parametrize the ALP 
mass as 

)- = tuq for T ^ A, 



mo[^] forT»A. 



/3 



(3.10) 



Here T is the temperature of the new sector. Naively it makes sense to assume A' ~ A" ~ A. 
At temperatures larger than A, electric-screening damps long range correlations in the plasma 
and thus the instantonic configurations, resulting in a decrease of the ALP mass. In specific 
models the exponent (5 can be obtained for instance from instanton calculations but here we 
will treat it as a free parameter. Assuming the onset of ALP coherent oscillations happens 
in the mass suppression regime, it is easy to obtain an expression for mo /mi which is the 
expected enhancement in the DM abundance. We find 

^ > (3 X 1.66^5^)2"+^ (3.11) 



(3.12) 




A" J 

and the factor that controls the enhancement is 



y^mompi A' /mpi 



Unfortunately these models can provide only a moderate enhancement of the DM density 
with respect to the constant m^ case. The gained regions for the ALP DM case for values 
of /? = 1,3,5,7,9 can be seen in Fig. 1 from bottom to top (the lowermost region mi = 
mo corresponds, of course, to /3 = 0, i.e. to the previously considered case). Actually, 
even considering unrealistic huge values of (3 does not help much, as can be seen from the 
asymptotic approach of the highest /? cases. This is refiected by the finite limit of Eq. (3.11) 
when /3 — )• oo, but it follows from its definition, Eq. (3.10). In the /3 — )• oo limit is a step 
function of temperature m^ = 0(T — A), and the relation mo = A^//^ determines A from 
mo and f^f, (or g) so each point in the m,^ — g parameter space has an implicit maximum DM 
abundance, independent of (3. 
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Surprisingly, it appears that the crucial assumption that leads to these conclusions is 
that A' ~ A because it does not allow to consider arbitrary small values for A for a given 
mass. Therefore, models in which A' ^ A imply generically higher DM abundance and 
therefore require smaller initial amplitudes (pi and consequently stronger interactions more 
prone to discovery. Unfortunately, at the moment we cannot provide for a fully motivated 
example. 

Finally let us note that the hidden sector responsible for the thermal mass of the ALP 
can have implications for cosmology if it survives until SM temperatures below the MeV 
range. For instance it can behave as dark radiation or dark matter during BBN. These 
constraints have to be studied case by case but can be easily avoided if the temperature of 
the hidden sector is smaller than the SM bath. 

3.4 Survival of the condensate 

The ALP CDM scenario can be tested via its coupling to photons. However, a blessing for 
detection purposes might also be the model's curse. Firstly, the two photon coupling endows 
ALPs with the decay channel 77. The corresponding lifetime in vacuum^ is 

r, = = « L3 X 10^^ s ( ,) ' . (3.13) 

ALPs with a lifetime shorter than the age of the universe, ~ 13.7 Gyr, cannot account for 
the DM observed in galaxies and have to be discarded. In the — g plane this corresponds 
to the region in the up-right corner of Fig. 1, labelled "r < lO^^s" and excluded from further 
discussion in this paper. Even if ALPs have much longer lifetimes, the few decay photons 
can still be significantly above the measured photon backgrounds. We explore this possibility 
later on. 

Secondly, ALPs from the condensate can be absorbed by a thermal photon which is 
either on-shell (7^ — )• 7*) or off-shell (7*^ — )• 7). Off-shell photons 7* are understood to 
be absorbed or emitted by another participating particle. The inverse-Primakoff process is 
a notable example of the latter, with the extra particle being a charged particle from the 
plasma, for instance an electron. 

The thermalization rate of ALPs due to the Primakoff process has been treated in the 
literature and found to be [46, 47] 

r<^,± = (log (^) + O.8I94) . (3.14) 

This is much faster than the decay rate ~ g'^m^, since at early times the temperature exceeds 
niif,. However, this rate is thermally averaged over ALP energies with a thermal distribution. 
If this rate would apply to the thermalization of the condensate, i.e. of the zero mode, the 
condensate would perish if g'^T^ ever exceeded the expansion rate H ~ T^/mpi. The higher 
temperatures would be the most relevant and the condensate survives if 5 < 1/-v/Tr"1pi- 
However, let us convince ourselves that this is not the case. 

First, note that the inverse-Primakoff process (/> + — )• 7 + is exponentially sup- 
pressed at high temperatures because the energy of the incoming electron E has to be enough 

^In the presence of a photon thermal bath the decay is stimulated by a factor 1/(1 — e"'"*''^^)^. 
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to produce a photon, which at finite temperature has a non-zero mass ~ eT. This thresh- 
old impUes 2m^E > or E/T > e^T/m^j) which can be huge if m,^ is tiny^. Since the 
abundance of these electrons will be exponentially suppressed so would be the rate. 

This is however not a showstopper, since the threshold can be easily overcome by con- 
sidering an additional photon in the initial state, i.e. 7 + 1^ + — 7>7 + e^. However, these 
and similar processes are also suppressed at high temperatures T ^> m^. The reason is that 
the two-photon coupling involves the derivative of the axion field (since F^^F^^^ is a total 
derivative) whose only component is the time component d^cj) ~ nrifffCj) (after ti) and therefore 
all amplitudes involving the zero mode ALP absorption are necessarily proportional to 
and absorption probabilities to m|. To provide a more complete example, we can compare 
the rate of Compton scattering of a photon of energy w, denoted by Tc{^), with the rate 
of the process ^{(jj) + i;^ + — t- 7 + e^. In the limit of very small the ALP absorption 
and the virtual photon scattering factorize and the differential ALP absorption rate due to 
photons of energy u is 

{2ujm^Y + {ujTc{uj)Y 

where dn^[uj) is the density of photons with energy uj and /3 ~ 1 is the photon velocity. We 
have included Tc as the imaginary part of the propagator to account for the finite photon 
lifetime in the thermal bath. For the rate of ALP absorption we have to integrate over the 
initial photon energies thermally distributed. In the — t- limit it is 

T,c-9'T'^ (3.17) 

where (Tc) denotes a certain thermal average, expected to be ~ a^T bearing phase space 
factors. In practice this is a very large suppression factor with respect to the naive thermal- 
ization rate g^T^, which leaves the parameter space shown in Fig. 1 untouched. We believe 
that a similar kind of reasoning can be applied to other ALP couplings of derivative type, like 
to fermions^. We plan to expand on these arguments at length in a forthcoming publication. 

If large scale magnetic fields are generated in the early universe, for instance during the 
electroweak phase transition, they can trigger ALP-photon oscillations. The ALP condensate 
then oscillates into a large scale coherent electric field which can be easily damped due 
to the high conductivity of the primordial plasma [49]. However, the high conductivity 
strongly suppresses the ALP-photon mixing and in practice the required magnetic fields are 
unreasonably large. For instance, if we assume magnetic fields are frozen into the medium, 
their strength redshifts as B = Bq{1 + zY = Bo^ao/aY and requiring no damping of the 
condensate produces the constraint 

(10-10 GeV-0 (3^) 109 GeV ^ ^ ^^'^^^ 



®One can check that after electron-positron annihilation, when the formulas we have used are not valid 
anymore, the evaporation is still slow. 

'^It seems that this argument was not taken into account in the axion case [48]. However, there it has no 
phenomenological consequences. 
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Figure 2. Region that would be affected by resonant ALP-photon oscillations leading to the evapo- 
ration of the ALP condensate in the case where primordial magnetic fields exist and they evolve with 
redshift as B = Bo{l + z)"^ . 



where Tg is the temperature of the universe when magnetic fields form. This damping would 
have implications in the ALP parameter space Fig. 1 if strong primordial magnetic fields of 
very early origin are eventually discovered and conversely, if WISPy ALP DM is established by 
direct detection experiments, one could use this bound to constrain the existence of magnetic 
fields in the very early universe. If the magnetic fields are produced during the electroweak 
phase transition, the plausible constrains are not very promising. One would then need to 
have primordial fields of the order of S ~ (close to equipartition with radiation), which 
means a value Bq ~ today, for them to have consequences for ALPs which are not 
already excluded by stellar evolution, i.e. for g < 10"^'' GeV~^. These fields are very close 
to the exclusion limits for primordial magnetic fields (if not excluded) and they would imply 
distortions of the CMB spectrum for sub-eV mass ALPs, see [50] and references therein. 

Note that the constraint (3.18) is independent of the ALP mass. This happens for 
relatively large values of Tb but if the primordial fields are created relatively late, a resonance 
in the ALP-photon oscillations (not considered in [49]) will dominate the ALP condensate 
evaporation, and this does depend on the ALP mass. This resonance is very sharp in time 
and only depends on the size of the magnetic fields at that particular time, and therefore is 
independent of the very early universe. We postpone the quantitative aspects to the next 
section, devoted to hidden photons, because HP-photon oscillations are the most relevant 
HP evaporation mechanism and the required formalism is essentially the same. The resulting 
bound is displayed in Fig. 2. 

3.5 Thermal population of ALPs 

The Primakoff process is able to create a thermal ALP population if r^^± exceeds the ex- 
pansion rate at some point in the history of the universe. This population competes with the 
condensate to form DM. The phenomenological implications of this population have been re- 
cently reviewed in [51] (see also [12]). In particular, thermal ALP CDM exceeds the measured 
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Figure 3. Exclusion bounds on axion-like particles from relic photons in the mass-lifetime parameter 
space. 

value Eq. (2.11) unless 

m<^<154ev(^°^), (3.19) 

where is the temperature at which r^^± ~ H where the ALP bath decouples from the SM 
bath. This bound assumes that 

r.~10^Gev('°"°°°^")' (3,20) 

is somewhat smaller than the reheating temperature. In the large mass region this puts a 
quite stringent limit on ALP CDM. However, it can be circumvented in scenarios with a low 
reheating scale and we therefore do not include it in Fig. 1. 

3.6 Detecting photons from ALP decay 

Even if the ALP lifetime is longer than the age of the universe some ALP decays inevitably 
happen and the resulting monochromatic photons can signal the existence of ALP CDM. We 
show the bounds coming from relic photon detection in Fig. 1 and more specifically in Fig. 3, 
where they are plotted in the — r parameter space (where r is the ALP lifetime). 

To obtain these bounds we have compared the extragalactic background light (EBL) 
spectrum from Ref. [52] with the estimated flux of photons due to the decay of the ALP [53, 
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54], 



dFE _ 1 p^{z)/m^ _ 

dEdn ~ 2-K H{z) (1 + zf 

. f^V'^xp f-^ i^T) e im, - 2Eo) (3.21) 

2tttHo \m^/2) ^\ t \m^/2) ] ^ ^ 

where we wrote all present quantities with a subscript "0" and the decay-redshift as 1 + z = 
{m^/2)/ Eq, with £"0 being the photon energy today. 

The EBL spectrum shows no particular features. Thus it can be used to exclude any 
ALP which would provide too strong signal. From the comparison we exclude the portion of 
parameter space labelled "EBL" in Figs. 1 and 3. 

Galaxies, being denser, should provide an enhanced signal for decaying ALPs. Again, if 
the decay rate of ALP CDM is high enough, we should be able to detect a spectral line whose 
energy is -E = m^/2. Axions in the visible part of the spectrum [55-57] and sterile neutrinos 
in X-ray spectra [58, 59] have already been searched using this technique. Following these 
references and conveniently rescaling their results for ALPs, we obtain the exclusion bounds 
in Figs. 1 and 3 labelled "Optical" and "X-Rays", respectively. The bounds from searches 
of gamma-ray lines in the data from FERMI [60] have been discussed in [51]. Unfortunately, 
the very suggestive yet tentative claim of a detection at Eq ~ 130 GeV [61], can only be 
attributed to annihilating DM and not to decaying DM. 

The search for ALP decay photons is blinded in the ultraviolet range by the humongous 
opacity of the atmosphere. However, a substantial amount of ionizing radiation at high 
redshifts can strongly affect the reionization history of the universe. The optical depth of 
reionization has been estimated by WMAP7 to be 0.088 it 0.015, out of which nearly a factor 
0.04 — 0.05 can be attributed to a fully ionized universe up to redshift ~ 6 as supported 
by the absence of Ly-a features in quasar spectra. The origin of the remaining part of the 
optical depth is still matter of debate, and thus we require the ionization caused by decaying 
ALPs to be less than what this fraction would require. Assuming one ionization per ALP 
decay photon, we have computed the reionizaton history using RECFAST [62] and obtained 
the modified optical depth from z = Q to z = 100, tq. In the region labelled "xion" in Figs. 1 
and 3, ALP decays produce too early reionization and exceed the measured tq ~ 0.04 by 1 
standard deviation. This is a conservative bound. If we assume that the full energy of the 
photon can be converted into ionization this bound strengthens increasingly with mass up to 
one order of magnitude at the largest masses for which ionization is effective, < 300 eV. 

3.7 Other indirect observational constraints on ALPs 

The propagation of photons and cosmic rays in today's ALP background might have poten- 
tially detectable consequences. The ALP background implies time-dependent birefringence 
which produces rotation of the polarization plane of photons propagating across the uni- 
verse [63]. For simplicity, let us consider the case in which the ALP field is homogenous and 
the frequency of light is bigger than the ALP mass, uj ^ m^. Then the rotation of the polar- 
ization plane in the WKB approximation is independent of uj and given by the line-of-sight 
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integral 



/\ip g I dl(pr^ ~~ (cos(m<^/ + /3) - cos(/3)) 

eV 



- 2 X 10~^ \^_,^^ , _i — (cos(m<^/ + /3) - cos(/3)) (3.22) 



where /3 is an unknown phase and we have neglected the expansion of the universe (including 
the expansion is straightforward). Clearly, this effect is 0{\) only for extremely small masses, 
and this only if we observe objects at distances ^ For instance, the bounds implied 

by the non-observation of such rotation in the emissions of AGNs [64] seem to be irrelevant. 
Indeed, recalling that todays ALP density, and therefore its field, is bounded by (2.12) and 
using (j)i = N9ia/2'Kg with A/" = 1 and Oi < tt one finds 



/ 1 8 X 10"^'^ eV 

A99< 10"%/^ (cos(m0^ + /3) -cos(/3)). (3.23) 

This shows that ALP DM can produce only minute polarization changes^. In (3.22) and 
(3.23) we have used the average value of DM which is much smaller than the values currently 
present in the galaxy or in the universe after recombination, which are bigger by factors ~ 10^ 
and 10^, respectively. Unfortunately, this does not change our conclusions. Even going to 
the largest possible distances, i.e. looking at the CMB, the effects have been found to be 
small in the range of masses relevant for ALP DM [7]. Other consequences, like the existence 
of gaps in the photon dispersion relation at frequencies close to integer multiples of rriff, have 
been discussed in the literature [66] and found to be not promising for the detection of ALPs. 

Another possible signature would be the radiation generated by cosmic rays propagating 
in the time-varying ALP background. Low energy radiation was considered in [67] and 
found to be many orders of magnitude below the backgrounds. At high energies, the inverse 
Primakoff process g + (/> — )> g + 7 could in principle contribute to the relatively small gamma- 
ray backgrounds, however the contribution arising from of inverse Compton scattering of 
cosmic rays off the CMB photons q + 7cmb —^9 + 7 dominates because of the much larger 
cross section. 



3.8 Direct experimental and observational constraints on ALPs 

So far we have focussed on the astrophysical and cosmological constraints specific to ALP 
CDM. However, in Fig. 1 we also show a variety of constraints that arise from the properties 
of ALPs alone without them needing to form all or part of the DM. In the mass region 
shown in Fig. 1 three constraints are most important. The constraint labelled "ALPS" is 
the result of the light-shining-through- walls experiment^ ALPS [69]. As we can see this 
experiment does not yet test the ALP CDM region. However, in the near future significant 
improvements by several orders of magnitude are expected by enhancing the signal with 

^This rotation is proportional to the amplitude of the sidebands looked for in the experiment proposed 
in [65] to search for meV mass axions or ALPs. Therefore, Eq. (3.23) seems to make their detection quite 
challenging. 

®In such an experiment strong magnetic fields are used to induce photon- ALP oscillations in incoming laser 
light. The very weakly coupled ALP state can then pass through a wall. On the other side the ALPs can 
then oscillate back into photons which can be detected. See [68] and references therein. For HPs the same 
technique works but no magnets are required. 
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resonant cavities [70, 71]. Moreover improvements are expected from upcoming experiments 
in the microwave regime [72-74] 

The "CAST+SUMICO" constraint arises from the hehoscopes^o [75] CAST [76, 77] and 
SUMICO [78]. As we can see these experiments already exclude sizable regions of the ALP 
CDM parameter space. Also for these experiments significant improvements are expected in 
the future. In particular if a next generation axion helioscope such as lAXO [79] is realized. 

Finally, the bound "HB" arises from comparing the observed cooling rate of horizontal 
branch stars with the expected rate. This places strong bounds on additional energy losses 
caused by a production of ALPs in the star's core [80, 81]. These bounds are currently 
the strongest and probe the ALP DM region. However, as these bounds are limited by 
astrophysical uncertainties we expect that the more controlled experiments discussed above 
will overtake them in the not too distant future. 

As already alluded none of these experiments make use of ALPs being DM. This makes 
them particularly model independent, but also ignores a potential plentiful source of ALPs. 
We will return to haloscopes which indeed exploit this source in Sect. 5. 

4 Hidden photons 

In a recent article [11], Nelson and Scholtz have considered the possibility that the misalign- 
ment mechanism could also be applied to generate a population of hidden photons (HPs), an 
Abelian gauge boson under which SM particles are uncharged. The Lagrangian is 

1 ml, 

/: = --X^,X^^ + ^X^X^ + £grav + /:/, (4.1) 

where is the HP gauge field and X^'^ its field strength. Moreover Ci contains the inter- 
actions with the Standard Model particles and £grav specifies potential non-minimal gravita- 
tional couplings discussed below. The HP mass might result from the Higgs or Stiickelberg 
mechanisms. In the first case, we have to worry when the phase transition happens and 
we might have a similar scenario to the one sketched in the previous section for Nambu- 
Goldstone bosons. Also a Higgs particle appears in the spectrum, with mass ~ yfXm,^! j 
where gh is the hidden sector gauge coupling and A the Higgs self-coupling. Even if we 
take gh to be really small, the Higgs particle phenomenology tightly constrains this scenario, 
especially for the sub-eV values of m^i we explore [82]. As in the original proposal [11], 
we focus therefore on the Stiickelberg case, which occurs naturally in large volume string 
compactifications [83-85]. In this case, there is no SSB phase transition. 

Let us briefly discuss the evolution of the HP in an expanding Universe^^. For reasons 
that will become clear at the end of this discussion we also include a non-minimal coupling 
to gravity of the form-*^^ 

^grav = '^RX^X^. (4.2) 

For simplicity let us focus on the homogeneous solution, diX^^ = 0. The equation of 
motion then enforces Xq = 0. As explained in [86] the invariant X^^X^ = —l/a?'{t)XiXi is a 

^"Helioscopes use the same idea as light-shining-through walls experiments. The ALPs are however produced 
inside the sun from photons interacting with the electromagnetic fields of electrons and ions in the plasma. 

We would like to thank Valery Rubakov and Christof Wetterich for noticing an error in the treatment of 
the cosmological evolution of a vector field (which is also present in Ref. [11]). 

^^We use a coordinates such that ds'^ — dt^ — a^(t)dxl, i.e. the metric is y^i/ — 
diag(l, —a?{t), —a^{t), —o?[t)). Moreover, the gravitational part of the Lagrangian is Cqyk. ~ — -R/(167rGjv). 
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coordinate independent measure for the size of the vector and it is convenient to introduce 
Xi = Xi/a{t). Using this the equation of motion is, 

Xi + mXi + (m^, + {l-n){H + 2H'^)^ X, = 0. (4.3) 

The energy density is 

p{t) = To=\ + m^yXiXi + (1 - K)H^XiXi + 2(1 - K)HXiXi^ . (4.4) 

For H <^ my and H <^ rriy the expressions (4.3) and (4.4) reduce to the same form as 
Eqs. (2.2) and (2.4), independent of the value of k. In consequence, the same approximate 
solution Eq. (2.10) holds for the HP case. In particular for my = const the energy density 
behaves (for any value of k) just like that of non-relativistic matter p{t) ~ l/a^(i). 

In the following we will derive constraints that arise from the evolution in this phase of 
the evolution. When the expansion rate is slow one can check that Xi actually is the properly 
normalized field. Therefore, from the following subsection on, we will simply drop the bar. 

Let us now have very brief a look at the initial conditions and the evolution before 
oscillations begin. In contrast to the ALP case, there is no natural value for the initial 
value of the amplitude^'^, Xi^i. Therefore choosing a suitable value for Xi^i any HP model 
can provide DM, i.e. nature might have tuned the value of Xi^i to fit the observed value 
Eq. (2.11). However, one might want to be a bit more ambitious and also consider the 
evolution before ti. Here, we will only consider the simplest 1. Then Eqs. (4.3) 

and (4.4) reduce to the same form as in the scalar case [86]. Accordingly, for H ^ m^' both 
the field Xi and the energy density are approximately constant, and we can think of these 
values being stuck there from some time at the beginning of, or even before, inflation. 

One further comment is in order. The direction of the X fleld remains unchanged for 
most of the universe history. However, it is conceivable that it changes during the process of 
structure formation where inhomogeneities in the gravitational potential grow relevant and 
Eq. (2.2) is no longer valid. This has important consequences for direct detection, that will 
be discussed in Sect. 5. 

4.1 Survival of the condensate 

Under the assumption that SM particles are uncharged under the hidden photon gauge group, 
the dominant interaction between the visible and the hidden sector is through gauge kinetic 
mixing between photons and hidden photons [87], 

1 1 rn^.,/ Y 

C = -^F^uF^'^ - -X^uX'^'^ + -^X^X'' - ^F^uX^'' + J^A^, (4.5) 

where Afj_ is the photon fleld and F^u the corresponding fleld strength, and is the current 
of electrically charged matter. Kinetic mixing is generated at one-loop by the exchange of 
heavy messengers that couple both to the ordinary photon as well as to the hidden photon, 
its natural value therefore being determined by the visible and hidden gauge couplings via 

^^If the hidden photon mass arises from a Higgs mechanism one may wonder if there is no hmitation on the 
field value from the fact that one component of the field arises from the "eaten" Goldstone boson which has a 
limited field range, similar to the situation we discussed for ALPs. Ifowever, in effect the extra field component 
of the gauge field corresponds to a derivative of the Goldstone boson, which is not bounded. Similar reasoning 
can be applied to the Stiickelberg case. 
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X ~ e(7/i/(167r^). In field theory and in compactifications of heterotic string theory, the 
hidden gauge couphng is of order one and thus x ~ 10^'^ [88, 89]. In large volume string 
compactifications, the hidden gauge coupling gh can be very small and there is no clear 
minimum for x- values in the 10~^^ — 10~^ range have been predicted in the literature [84, 85]. 

For a Stiickelberg mass the same string scenarios typically prefer values > 10~^ eV. 
Nevertheless, to be as inclusive and model- independent as possible we will take both x and 
niY as free parameters in our phenomenological study. 

By means of the re-definition = A^ — xX^, = X^ we can identify the propagation 
basis in vacuum where the kinetic mixing has been removed and the Lagrangian 



The universe is not empty though. The interactions of photons with the charged parti- 
cles in the primordial plasma induce refraction and absorption. This can be described with 
an effective photon mass squared = + iuV, where both the plasma mass and 
the interaction rate T only depend^'^ on the photon frequency, w, and the modulus of the 
wavenumber, k [90]. The corresponding effective term in the Lagrangian can be written as 
C = M"^ A^A^ /2. This is non diagonal in the {^1, X} basis and crucially determines the 
propagation eigenstates and their eigenvalues [90-93]. 

grows with the charge density nq and the temperature and therefore varies signif- 
icantly during the history of the Universe. After inflation we can assume the universe has 
been efficiently emptied and thus = 0. During reheating, uq and T grow very rapidly 
up to extremely large values and then decrease with smooth power laws until today, where 
we shall assume them to be zero again. Forgetting pre-/reheating times for a moment, at 
sufficiently early times is the largest scale in the problem. Neglecting H, the e.o.m. for 
the photon A = A — xX is that of a damped oscillator, while the orthogonal combination 
X + x-^ is (almost) a free field with a mass ~ m^'. The latter is the combination which can 
form a condensate. 

Today tends to 0, so A and X are the propagation eigenstates, their e.o.m. decouple 
and both follow the solution outlined in Sect. 2. Therefore, the condensate has to follow a 
trajectory in field space from X + x^ to X without losing too much amplitude. 

Assuming that the universe expands slowly compared to the microscopic timescales of 
photon absorption {dV/dt <C F^, or F S> -ff) the two eigenstates adiabatically follow their 
static solutions. The adiabatic propagation eigenstates are given by an effective mixing angle, 
Xeff, the angle between the sterile state {X + xA) and the condensate. The eigenstates decay 
with rates Fi = (1 — Xeff)^ and F2 = X^q^- At the lowest order in x, the effective mixing 
angle can be approximated by 



where = maxjxmy , m^'F}. 

Before recombination is positive so there might be a moment where the mixing is 
resonant. Characterizing it by the temperature of the universe at the time of the resonance 
m^(rres) = my we can distinguish three regimes 

^*In an isotropic homogeneous universe. 



looks like 



C = --F,,,F^'^ - -X^uX^^- + -^X,X^ + [a^ - xX,) . 



(4.6) 



Xeff - 




(4.7) 



- 20 - 



• At high temperatures T ^ Tj-^s-' rn-y ^ my and the mixing is suppressed as XeS — 

• At low temperatures T <^ T^cs^ m-y <^ my and we recover the vacuum mixing parameter 
Xcff — X (since usually < my). 

• At resonance T = T^es- the mixing is significantly enhanced. If X"i7 > T the mixing 
angle is maximal and one can have a resonant transition similar to the MSW effect in 
neutrinos. For small x the mixing still is enhanced with respect to the vacuum case by 
a factor m^/T ^ 1. 

Thus, in general the condensate does not move smoothly in field space from X + x^ to 
X. It overshoots X during the resonance only to then return to X. 

Let us first discuss the case of over-damped oscillations during the resonance, i.e. Xeff ~ 
X^w,y/r^. The amplitude of the HP condensate decreases by the ratio 

^ = {-^)"\^ f 1 ,,r.) . (4.8) 

^2,initial \ C^today / \ ^ ^tinitiai / 

Approximating T/Tq = ao/a, we can write T2dt as T2/Hd\ogT. The function T2/Hx^ = 
rXgff/-ffx^ is plotted for different HP masses in Fig. 4. In the plot we have used 

M2 = luja = ujI vT + ^^^Th 4.9 

1 + [ujtY 

where u = my is understood, wp is the plasma frequency, Fxh = o'Th'^e is the absorption 
due to Thomson scattering off the electrons and r the characteristic time between electron 
collisions, which sets the friction time scale^^. 

As expected the decay rate is heavily enhanced at T^es-, and rapidly drops at higher 
temperatures. The integral can be well approximated by the contribution near the resonance. 
Using m^ = my + |(imy /dT|i.es(T' — Tres) and evaluating the other T-dependent quantities 
at Tres we have 

1 /■ "^initial dTT2 2 ^ "^7' _ ''"2 



2 Jr , T H ^ ^ 2 rH,^^ 2 ^^'^^^ 

•'-'today 

where r = dlogm^/dlogT is an 0(1) factor, see Fig. 2 of [90]. Moreover, exp(— T2) gives the 
damping of the total energy density due to the resonance. Note that the result is independent 
of r. 

If the resonance is not over-damped, i.e. x^y > T, the HP survival probability after 
the level-crossing can be approximated by the Landau-Zener expression used in [93]. The 
result coincides almost exactly with (4.8), only requiring an extra factor of 2 in the argument 
of the exponential, i.e. T2 — )• 2r2. 

The plasma effects were neglected in Ref. [11] which forced the authors to conclude that 
the most relevant period for the evaporation of the condensate was at high temperatures 
(the corresponding "approximation" is shown as a dotted line in Fig. 4). As we have seen, 



^^The zero-frequency conductivity ao = ujpT is normally used instead of r. We have followed the approx- 
imations for (To that can be found in [49] and supplemented them in the non-relativistic electron case with 

2 

Thomson friction ao ~ ^^f^ " where a is the fine-structure constant, and ne,me the electron density and 



- 21 - 




T[keV] 

Figure 4. The decay rate of the HP condensate, normahzed to Hx^ , as a function of temperature, 
for different HP masses. The curves, ordered by their resonance from left to right, correspond to 
m^/ = 10~^, 10~^, 10"'^, 10~^ and 0.1 eV. The dashed hne corresponds to the expression used in 
Ref. [11]. 



at high temperatures the evaporation process is strongly suppressed by the effective mixing 
angle. The resonance dominates the condensate evaporation and for the sub-eV masses of 
main interest here happens at temperatures smaller than the electron mass [92]. 

Before we continue, it is now time to describe the differences between HP-photon and 
ALP-photon oscillations in the presence of a background homogeneous magnetic field of 
primordial origin. The physics is essentially the same, except that for ALP-photon oscillations 
the mass mixing term in (4.7) is gBm^f, [49] with B the magnetic field strength, instead of 
X'm^i'^. The important difference is that the value B is expected to change during the 
evolution of the universe. The most relevant effects appear when magnetic fields are frozen 
to the medium, in which case B = BQ(ao/a)'^ ~ Bq(T /Tq)"^ . In this case, the integral (4.10) 
tends to be dominated by the highest temperatures, which in this case are the temperatures 
at which the magnetic fields are created Tg 

/ \ 2 / n \ 2 

r2 



2 



ALP VlO-i°GeV~V V3nGy IQSGeV ^ ' ' 



Imposing T2 < 1 leads to the bound (3.18). The resonance contribution can be computed, in 
full analogy to (4.10). We obtain 



T2 

2 



ALP 2 rm^Hj 
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Let us now come back to the discussion on HPs. Since we can adjust the value of 
,initiai almost without restrictions^^ it seems that we can overcome any evaporation factor, 
even if it is enormously exponentially suppressed, and still have today HP CDM. There are, 
however, several limitations to this due to the fact that the evaporation process dumps some 
energy into the photon bath, 



This photon injection dilutes neutrinos and baryons with respect to photons, which can be 
constrained by the effective number of relativistic neutrino species A^^''^ and by the successful 
agreement of CMB and big bang nucleonsynthesis (BBN) to estimate the baryon to photon 
ratio r]B- Considering three massless neutrino species, a number of effective neutrinos smaller 
than 2.39 is excluded at 95% C.L. [94], which translates into a photon temperature increase 
T'/T ~ 1.06. A similar value is obtained from BBN limits but the corresponding bound 
only applies to injections between CMB and BBN times, while the N^^ one applies until the 
neutrino decoupling T ~ 2 MeV. If the injected photons thermalize, the photon temperature 
increases to 



where we used pcBM ~ ^mpriBn^, with nip being the proton mass and n.y the photon number 
density. We see from Eq. (4.14) that to obtain a significant increase in temperature one needs 
T2 ^ — ln{r]B) ~ 21. Imposing^^ T' /T^^s < 1-06, we can exclude the region above the curve 
labelled ''Nf" in Fig. 5. 

If the resonance happens below a critical temperature, the interactions of photons with 
the relic electrons and ions of the plasma might not be enough to fully recover a blackbody 
spectrum. The photons are initially injected at ultra-low energies, oj ~ wpi <^ T, at which 
inverse bremsstrahlung (7 + e~ + — )• e~ + p'^) and, to a lesser extent, inverse double 
Compton scattering (7 + 7 + 6" — )• 7 + e~) are always efficient in absorbing them. This 
results in the establishment of a blackbody distribution at a higher temperature than the 
initial, but only at low frequencies, where these processes are effective. The equilibration of 
the low energy and high energy parts of the spectrum can be achieved by a combination of 
Compton scattering (which is efficient in reshuffling photons up and down in energy) and 
the photon- number changing processes mentioned above which adjust the number density 
towards a blackbody. The injection of a relatively high number of low energy photons, 
5n^/n^ ^ 5p^/p^, has not been considered in the literature and the use of numerical methods, 
even if possible, is clearly beyond the scope of this paper. In order to obtain a first educated 
guess, we have followed the analytical solutions derived in the limit of small distortions 
in [95] and the powerful constraints set by FIRAS on a possible chemical potential /x and a 
Comptonization y distortion which imply bounds up to Ap^/p^ < 0(10"^) [96], depending 
on Tres- These arguments lead to the constraint shown in Fig. 5 labelled "CMB distortions", 
which is the strongest requirement in the HP mass range neV — 2 x 10~^ eV. 

The lower limit of the CMB constraints, my ~ neV, corresponds to resonances hap- 
pening around the onset of recombination. Smaller HP masses suffer resonant transitions 

^^One may wish to impose that the initial field value is smaller than mp\. But this does not lead to any 
relevant additional constraints. 

^^In this high T2 regime, our formula (4.10) is not completely consistent. But, as long as T'/T ~ 1, it should 
provide a reasonable estimate. 



Ap = /OCDM(e'^^ - !)• 



(4.13) 




(4.14) 
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around this epoch because neutral Hydrogen contributes to with a negative quantity 
which increases with time and eventually makes = by compensating the contribution 
of free electrons. HPs with masses down to ~ 2 x 10~^^ eV have their resonance around this 
epoch [93]. It is very likely that CMB distortions will be produced for sub-neV masses, but 
since we cannot ascertain the fate of the low energy photons injected during recombination, 
we cannot use the results of [95] and we must leave the study of this region for a future 
detailed numerical analysis. It is intriguing that a possible small distortion of the CMB spec- 
trum due to HP CDM can in principle appear in the next generation of CMB spectral probes, 
such as PIXIE [97]. We believe, these signatures can be quite distinctive and strengthen the 
physics case for such probes, since they can help identifying the nature of the DM. 

There is however a simple way of constraining the kinetic mixing of CDM HPs in this 
low energy region. At late recombination, the density fluctuations are already imprinted in 
the CMB and they provide us with the estimate (2.11) for the DM density. Since this value 
agrees roughly with the average DM density observed today we cannot allow the resonant 
depletion of HPs into photons at this epoch or any later, thus requiring r2 < 1. This bound 
is depicted also in Fig. 5 and labelled as "t2 > 1". In order to compute it, we have used the 
model for ni^ provided in [93]. 

Hidden photons with masses below twice the electron mass, can only decay via an 
electron loop into three photons. This decay rate is extremely suppressed for low mass HPs, 

oc a^m^^, /m^ [90]. However, for the most massive HPs considered here the decay can 
be effective. Imposing that the population of decay photons is smaller than the diffuse X- 
ray backgrounds one can constrain the cold HP population. This was done in [90] and we 
reproduce their bounds, labelled "X-rays" , as the rightmost boundary of our HP CDM region 
in Fig. 5. 

4.2 Thermal population of HPs 

During the resonance a thermal population of HPs can be generated without interfering with 
the arguments above. The role of this population as DM has been addressed in [90] (which 
excludes the region labelled "Cosmology of thermal HP DM" in Fig. 5) and, as dark radiation, 
i.e. contributing to the number of effective neutrinos, in [92]. This conversion also produces 
distortions of the CMB spectrum, which can be constrained by FIRAS. These arguments 
give rise to the bounds labelled "CMB" in Fig. 5. 

These thermal contributions are very much constrained by astrophysics and cosmology 
and on the verge of detectability by solar HP searches or searches for the contribution of the 
7' — 37 decay to the interstellar diffuse photon background. 

One could ask the question whether it is possible to have HP CDM and also explain 
the small excess of A'^^^ observed by WMAP7 and other CMB and large scale structure 
probes. The situation is a bit tricky since, for the parameters required to have the latter 
effect (m-y' ~ ©(meV), x ~ 10~^) [92], the CDM evaporation is huge, T2 ^ 20, and the 
required absorption of the CDM energy decreases N^^ opposing the first effect. Further 
studies are required to see if such a scenario is plausible, but it seems complicated to avoid 
the BBN bound coming from the dilution of the baryon density. 

4.3 Indirect observational constraints on HPs 

For ALPs we have found powerful constraints from ALP decay and one may wonder if similar 
constraints exist for HPs. Indeed, as mentioned above, HPs can decay via 7' — )• 37 and this 
constrains the high mass region in Fig. 5. 
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Figure 5. Allowed parameter space for hidden photon cold dark matter (HP CDM) (for details see 
text). The exclusion regions labelled "Coulomb", "CMB", "ALPS", "CAST" and "Solar Lifetime" 
arise from experiments and astrophysical observations that do not require HP dark matter (for a 
review see [38]). We also show constraints on the "cosmology of a thermal HP DM". Note that 
only constraints on HPs with masses below twice the electron mass are shown since otherwise the 
cosmological stability condition requires unreasonably small values of the kinetic mixing, x- The 
four constraints that bound the allowed region from above, "t2 >1", "CMB distortions", "A^^^" and 
"X-rays" are described in the text. 

Similar to the ALP case one may also wonder about constraints from photon and cos- 
mic ray propagation. Photon propagation is essentially unaffected by a HP dark matter 
background since the combined photon-HP equations of motion are still linear and conse- 
quently the superposition principle holds. In other words photons pass right through the HP 
background without interacting. 

Cosmic rays on the other hand could scatter of the HPs via Compton scattering, q+j' — t- 
(7 + 7. However, for the relevant values of the kinetic mixing parameter the cross section is 
too small to have a significant effect. 

4.4 Direct experimental and observational constraints on HPs 

In Fig. 5 we have also displayed the existing experimental bounds on the existence of HPs 
which do not rely on HPs being DM. The bounds labelled "Solar lifetime" and "CAST", 
coming from the non-observation of HP emission from the Sun, exclude a large portion of 
parameter space [91]. It is clear that improving the sensitivity of future searches of solar HPs 
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one has access to new parameter space in which HPs can be CDM [98]. The solar hidden 
photon search (SHIPS) [99] in the Hamburg Observatory is currently exploring the sub-eV 
mass region greatly improving over the previous CAST experiment and will soon publish 
results. 

Light-shining-through-walls experiments are also a powerful tool in the search for hid- 
den photons [100], the current best bound being provided by ALPS [69], and shown in Fig. 5. 
As in the case of ALPs they currently do not probe the DM region, but significant future 
improvements are expected [68]. Indeed, here the microwave regime is particularly promis- 
ing [73, 74], with several new experiments already taking data [101, 102]. 

Finally, tests of Coulomb's law provide strong constraints in the low mass region, see 
Fig. 5. Although these bounds do not reach into the HP DM region one can hope that 
improvements can be made since the best experiments [103] on this are more than 40 years 
old. 



5 Direct searches with haloscopes 

One well-known tool to search for axion dark matter are so-called axion haloscopes [75]. 
Let us briefly recap the basic principle of a haloscope. Using the abundance of axions all 
around us the task at hand is to exploit their coupling to photons and to convert those axions 
into photons which can be detected. For axions this can be achieved by utilizing off-shell 
photons in the form of a strong magnetic field. Moreover, this conversion can be made more 
efficient by employing a resonator, resonant at the frequency corresponding to the energy 
of the produced photons. The energy of the outgoing photons is equal to the energy of 
the incoming axions. As the DM particles are very cold their energy is dominated by their 
mass. For axions in the natural dark matter window this mass is in the 1 — 100 fieV range 
corresponding to microwave photons. A number of experiments of this type have already 
been done [104-108] and further improvements are underway [109-111]. 

5.1 Axion-like particles 

For axion-like particles the experiment proceeds exactly as in the axion case [75] (and very 
similar to the hidden photon case described below). The important point is, however, that 
now even bounds which do not reach the predictions of axion models (axion band in Fig. 1) 
become meaningful since they test viable models. 

For completeness we recap. In the axion case, the power output of a cavity of volume 
V, quality Q and coupling k to the detector is 

-Pout = Kg'^V\Bo\'^ poGs,xion — Q , (5.1) 

with pq the local axion CDM energy density, |Bo| the strength of the magnetic field and 



|Bo|2y / dy|Ecav|2' 



For cylindrical cavities in the TMqio mode as, e.g., used in the ADMX experiment, ^/axion = 
0.68. The currently excluded region from various axion haloscopes is shown as gray area in 
Fig. 1. It already excludes a part of the ALP CDM parameter space. 
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5.2 Hidden photons 

Microwave cavity experiments looking for relic axions could also be used to constrain and 
search for the hypothetical cold HP condensate that we have discussed above. Starting from 
Eq. (4.5) we can follow the usual route and trade the kinetic mixing term for a mass mixing 
by performing a shift of the HP field X ^ X — x^. Neglecting terms of second order in x 
the equation of motion for the photon field A then reads, 

df^d'^A-' = xm^yX". (5.3) 

We can therefore see that the hidden photon field acts as a source for the ordinary photon. 

Let us first determine the strength of this source. As discussed above we can take 
X^ = 0. For the spatial components we write X. The energy density in the hidden photons 
is equal to the dark matter density. Therefore we have, 

2 

in 

PCBM = (5.4) 

At this point a comment concerning the direction of the HP field is in order. In the 
discussion above it was assumed that the HP field is homogeneous in space and points in 
the same direction everywhere. Due to structure formation the DM is clumped. A common 
estimate for the DM energy density on Earth is therefore the energy density in the galactic 
halo, 

0.3 GeV keV , , 

PCDM.local ~ n — > PCDM,average ~ n- (5.5) 

cm'^ cm'^ 
As we can see the local density is much higher than the average density, signifying that 
structure formation clearly is very important. This raises the question whether structure 
formation also infiuences the direction of the HP condensate. To answer this question one 
would have to study structure formation for a vector field like our HP field. This is beyond 
the scope of the current paper. Instead we will consider two possible scenarios: 

(a) The direction of the HP field is (essentially) unaffected by structure formation and all 
HPs point in the same direction (at least for a suitably big region of space). 

(b) The direction of HPs behave like a gas of particles with the vector pointing in random 
directions. 

In scenario (a) the HP direction is characterized by a fixed vector n, whereas in the case 
(b) we have to average the final result over all directions for h. With this understood, let us 
write, 

X(x) = n^^, (5.6) 

my 

with pq the dark matter density on earth. 

Let us now return to our cavity experiment. The photon field A inside the cavity can 
be expanded in terms of the cavity modes, 

A(x) = J]aiAr(x), I d3x|Ar(x)|2 = C„ (5.7) 
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with Ci the normaUsation coefficients. Using this expansion and inchiding losses in the cavity 
we obtain for the expansion coefficients, 



^ + + ^0 I "hO = Oiexp(-iwt), (5.8) 

with loq the frequency of the cavity and Q its quaUty factor. The driving force bi can be 
written as 



dVA*{x) ■ X(x) (5.9) 



and the frequency is given by the energy of the HPs, 

to = Ey my. (5.10) 
The asymptotic solution for the cavity coefficients is then, 

ai{t) = Uifi exp(-iwt) = j — exp(-ia;t). (5.11) 

Finally, the power emission of the cavity is related to the energy stored and the quality 
factor of the cavity as 

U 

Vout = i^Q^o, (5-12) 



where k is the coupling of the cavity to the detector and 

J ^i3x|Ar(x)p. (5.13) 



2 

Replacing in this equation and evaluating at resonance ujq = fn-y' we find 

Vont = t^X^m^^pQVG, (5.14) 

where the geometric factor Q is defined as 

|/dyA*'=^^(x) -nl^ 
^ ~ y /d3x|A--(x)|2 • ^^-^^^ 

In a cavity E = a; A. This geometry factor has exactly the same form as in the axion 
case (cf. [107] and the subsection above) but with the direction n replacing the direction of 
the external magnetic field Bq in the axion case. Accordingly, 

^ = ^axionCOS^(6'), (5.16) 



^^As in the axion case using this condition requires that the Q is not too large since after structure forma- 
tion the DM particles move with different velocities of order 300 km/s~ 10~^c. This restricts the maximal 
usable g to be < 10^ If axions form a Bose-Einstein condensate (see [16, 17] and Sect. 2) or the galactic 
exhibits special structures as suggested in [112, 113], the velocity distribution could have significantly nar- 
rower structures potentially allowing to utilize a cavity with higher Q. Moreover, in this case one could also 
benefit from using narrow bandwidth techniques. Strong magnetic fields typically limit the maximal Q < lO''. 
Therefore, the former only works for hidden photons, where no magnetic field is needed and one can use super- 
conducting cavities. The latter however, would work for both HPs and ALPs and has indeed been employed 
in ADMX [114]. 
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where 6 is the angle between the magnetic field direction used in the axion case (this is 
usually chosen to be optimal) and the direction n of the hidden photon field. 

We can now use this formula to constrain the HP CDM with the present microwave 
cavity searches for axions [104-108]. 

For scenario (a) the result depends on the relative orientation of the cavity at the time 
of the measurement with respect to an a priori unknown direction of the HP field. This 
requires detailed knowledge of the timing of the experiment etc. A conservative estimate for 
the sensitivity can be obtained, however, by assuming that all directions in space are equally 
likely and taking a value for cos^(0) such that the real value is bigger with 95% probability. 
The corresponding value is cos^{6) = 0.0025. For situation (b) one can average over all 
possible directions and obtain (cos^(9)) = 1/3. 

The results of the analysis are shown in Fig. 5, under the label "Haloscope Searches". 
Blue shows scenario (a) with our simplified estimate and red scenario (b). With a more careful 
analysis of the experiment, including the timing information, the sensitivity in scenario (a) 
would be very similar to the one in scenario (b). 

5.3 Non- resonant searches 

Both for ALP DM and for HP DM the coupling to photons can be orders of magnitude bigger 
than in the axion case. At the same time it is desirable to explore a bigger mass range. For a 
first broadband search one may therefore be prepared to give up the resonant enhancement 
in the cavity. Indeed the formulas in Sect. (5.2) up to Eq. (5.13) are valid also off resonance. 
One can obtain the general result by replacing 

In principle one can therefore explore a wide range off masses with a single fixed frequency 
cavity (at the price of some loss of sensitivity). However, it should be noted that in order to 
do so one also has to have a receiver sensitive in a sufficiently wide range of frequencies and 
also check for signals off the resonance of the cavity. 

6 Conclusions 

Vacuum misalignment in the very early universe is a very generic mechanism capable of 
producing a cold dark matter condensate from any very light, but massive (quasi-)stable 
bosonic field. Clearly, its self-interactions and its interactions with the SM should be very 
weak if it is to survive and to be a candidate for the observed cold dark matter today - 
in other words their particle excitations should be very weakly interacting slim particles 
(WISPs). 

We found that for the most prominent and theoretically well-motivated WISP candi- 
dates - axion-like particles (ALPs) and hidden photons (HPs), whose dominant interaction 
with the standard model arise from couplings to photons - a huge region in parameter space 
spanned by their masses and their photon coupling can give rise to the observed cold dark 
matter. 

WISPy CDM coupled to photons is particularly interesting because it can be probed 
by both astrophysical observations as well as laboratory experiments. Figs. 1 and 5 show 

^^Note Q may be a function of u. 



(5.17) 
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that sizable regions of this parameter space have already been excluded by experiments 
searching for WISPs without relying on them being dark matter. In particular, helioscopes — 
experiments looking for WISPs produced in the sun - have already probed regions that allow 
WISPy DM and will improve further in the near future. Soon also purely laboratory-based 
light-shining-through-walls experiments (in the optical as well as in the microwave regime) 
will reach a sensitivity that will allow them to test hidden photon DM with masses in the 
/LteV — meV region. 

The possibility that WISPs can form all or part of DM allows for additional search 
strategies. One way to exploit this is via haloscopes which already test both ALPs and 
hidden photons extremely sensitively in the //eV mass region. Moreover, searches for light 
emitted from decaying ALPs already provide interesting constraints in a large range of masses 
from mcV to MeV and will hopefully further improve in the future. 

Despite these bright prospects for future searches it is clear that huge areas of parameter 
space allowing for WISPy CDM are so far unexplored and demand not only new searches 
but also new search strategies - a challenge for both experiment and phenomenology. 
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